Strongly interacting bosons in a ID optical lattice at incommensurate densities 
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We investigate quantum phase transitions occurring in a system of strongly interacting ultracold 
bosons in a ID optical lattice. After discussing the commensurate-incommensurate transition, we 
focus on the phases appearing at incommensurate filling. We find a rich phase diagram, with 
superfluid, supersolid and solid (kink-lattice) phases. Supersolids generally appear in theoretical 
studies of systems with long-range interactions; our results break this paradigm and show that they 
may also emerge in models including only short-range (contact) interactions, provided that quantum 
fluctuations are properly taken into account. 

PACS numbers: 



The rapid progress in trapping and cooling atoms has 
rendered the study of "tailor-made" low-dimensional (D) 
systems [T] experimentally accessible. Both the dimen- 
sionality and the interactions can be controlled, allow- 
ing great flexibility in realizing almost arbitrary strongly- 
correlated physical systems. A superfluid-Mott insulator 
(SF-MI) quantum phase transition, driven by increasing 
the potential depth of the optical lattice (and hence the 
relative strength of interactions) beyond a critical value, 
has been observed for bosons loaded into an optical lat- 
tice in 3D 0, 2D [3], and ID [J. In addition, the Tonks- 
Girardeau gas, where bosons avoid spatial overlap and 
acquire fermionic properties due to strong repulsive in- 
teractions, has been experimentally realized in ID [5]. 

Recently, a new type of quantum phase transition was 
observed in ID in the very strongly interacting regime: 
for an arbitrarily weak optical lattice potential commen- 
surate with the atomic density of the Bose gas, a quan- 
tum phase transition into an insulating, gapped state, 
was observed, with the atoms pinned at the lattice min- 
ima H]. Theoretical studies of ID systems based on the 
sine-Gordon model indeed predict that above a critical 
interaction strength, the SF should become a MI even 
for a vanishingly weak optical lattice • 

In this Letter, we show that another interesting regime 
can be reached if the density is incommensurate with 
the optical lattice. The system is then described by a 
driven sine-Gordon model. In this model, the appear- 
ance of superfluidity (off-diagonal correlations) may be 
driven in two different ways, either by tuning the inter- 
action strength at constant lattice depth and commensu- 
rate period, as already realized experimentally [5], or by 
tuning the density or lattice parameter away from com- 
mensurability. We study the excitations of the system 
in the incommensurate phase and show that a supersolid 
(SS) phase may arise. In addition, for sufficiently large 
lattice strengths, a solid (S) phase is stabilized even at 




FIG. 1: T — phase diagram for a filling factor slightly 
above unit filling. SF, SS, and S indicate superfluid, super- 
solid and solid phases, respectively. The inset shows the same 
phase diagram but in terms of 7, the dimensionless interac- 
tion strength, instead of the Tomonaga-Luttinger parameter 
K. 

non-unit filling. 

Our approach to studying this system is as follows: 
First, following previous work [S{8|, we formulate the 
underlying ID interacting boson problem in terms of 
a quantum sine-Gordon field theory, with the devia- 
tion of the number density from commensurate values 
driving the appearance of kinks in the field. We care- 
fully obtain the density threshold for the commensurate- 
incommensurate MI-SF transition, finding that the SF 
appears for arbitrarily small deviations of the den- 
sity from the commensurate density, in agreement with 
Ref. 0. 

We next study the incommensurate regime, in which 
the excess particles appear as kinks of the sine-Gordon 
field, having an effective mass and effective interactions 
different from the bare particles. We extract these two 
parameters exactly from the underlying field theory and 
finally apply a functional renormalization group (RG) 



2 



approach to the path integral formulation of the many- 
body statistical density matrix to obtain the ground state 
properties of the system. The RG transformation shows 
that quantum fluctuations renormalize the interactions 
between the kinks to a power law form; this maps the 
system to the Calogero-Sutherland model 10 , an exactly 
solvable model of ID particles with long-range interac- 
tions. This finally allows us to propose a phase diagram 
for the incommensurate regime (see Fig. [T]). When the 
lattice potential is strong enough, the system solidifies. 
However, the S phase here is actually a lattice of kinks, 
and the number of particles per site is not fixed. At 
intermediate values of the lattice potential, we predict 
the emergence of a SS phase. SS phases usually occur 
in model Hamiltonians which include long-range inter- 
actions, and have a characteristic wavelength which is 
an integer multiple of the lattice spacing [TI]. The most 
striking feature of the system studied in this Letter is that 
while the original Hamiltonian contains only local interac- 
tions, the SS phase appears due to the finite-range nature 
of the interaction between the excitations. In addition, 
the periodicity of the SS phase found here is unrelated to 
that of the lattice, a qualitatively different behavior from 
the situations usually found in the literature. 

The microscopic description of a trapped gas of cold 
bosons in ID with contact interactions and in the pres- 
ence of a single-particle potential V{x) is 



H = 



dx 



^P\x) + V{x)^ 



dx -0^ {x)4'{x)ip{x) 



(1) 



where g is the strength of the (5-function interaction, ?/; 
{ip^) are bosonic annihilation (creation) operators, and 
m is the atomic mass. The parameter characterizing the 
strength of the interactions is the Lieb-Liniger parameter 
7 = mg/h^no, where hq is the average density. 

Writing = \J n{x) exp{—i9{x)), with n{x) the 

density and 9{x) the (real) phase, and using the Pois- 
son summation formula, the density operator may then 
be expressed as [U [8] 



n{x) = (no - -dx<j){x)j ^ 



^2ip[7TnQX-4>{x)] 



(2) 



where / dxdx4>{x) = 0. Eq. (2) yields an expression for 
the bosonic operators in Eq. ( 1 ) in terms of the new field 
(p. The appropriate bosonic commutation relations are 
satisfied if [dx4>{x),9{x')] = —iTrh6{x — x'); that is, 9 and 
dx4>/T^ are canonically conjugate variables. From Eq. ([2]), 
it follows that kinks in the 4> field correspond to particle- 
like excitations. This fact will be of great importance to 
us later on. 

In the long- wavelength limit, and in the presence of an 
optical lattice creating a single-particle potential V{x) = 




FIG. 2: Q > 0; the dashed line is dx<l>/T^, the solid line is the 
periodic potential, and the corresponding particle positions 
are indicated by dots. Notice how kinks (indicated by local- 
ized deviations from a straight line for 8x4'/''^) correspond to 
particles in excess of the commensurate particle density. 



(y/2) cos(47ra::/AL), the system of Eq. ([!]) may be de- 
scribed by an action of the form 



S[4>] 



dr / dx 

i-oo 8nK 



dr 



dx cos [4>{x) — Qx] , 



(3) 



where we have now seth= 1, scaled lengths such that the 
speed of sound is unity, and finally scaled (j) ^ (j)/2. Here, 
j3 = l/ksT, u = n^V , while K is the usual Luttinger 
liquid parameter. For bosons interacting via contact po- 
tentials, K may be expressed in terms of 7; for large 7, 
K (l-t-2/7)^, while for smaller interaction strengths 
7 it is given by tt/ — 7^/2/(27r). Notice that 

K 1, as it should for bosons with local interactions. We 
have also only kept the most relevant (least quickly os- 
cillating in space) terms and written Q = 27r (no — 2/Al) 
as the deviation of the average density from its commen- 
surate value. 

In the zero-temperature, /? — ^ 00 limit, Eq. ([s]) is for- 
mally equivalent to the model studied in Ref. [12 . It is 
also related to previous work on quantum Hall bilayer 
systems [131 HI] : with the important difference that the 
boundary conditions in the present case are / dx = 0, 
while in Refs. [T^[T5] (amongst numerous others), there 
is no such restriction on (p. This is crucial to the position 
of the boundary of the commensurate-incommensurate 
transition, and is due to the fact that we are working at 
fixed particle number [9]. 

Since kinks correspond to excess particles above the 
commensurate density (see Fig. [2]), fixing the particle 
density must fix the number density of kinks uniquely. 
But from its definition, Q is directly proportional to this 
excess particle density, so that the kink density must be 
proportional to Q itself. This immediately implies that 
Qc = 0, at least at zero temperature ^2E' . Mathemati- 
cally, this is a consequence of the boundary condition at 
the edges of the system, which implies that the commen- 
surate phase cannot exist unless Q = 0. For any Q > 0, 
a finite density of bosonic kinks appears. 
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The full analytical expression for the bare interaction 
between two kinks in a sine-Gordon model is given in 
Ref. [16j; its limiting forms are 



Vbarc(?') -5> 



^ exp(- 
V(2r), 



r > ^ 



(4) 



The effective width of the kinks is ^ = Xj \/2'KKnQV , so 
that this potential amounts to an impenetrable core plus 
a finite-range repulsion. 

The effective kink (and antikink) mass a is strongly 
renormalized by quantum fluctuations |17| . It may be 
obtained using a flow-equation RG scheme as described 
in Ref. 18; the salient points are that it is proportional to 
u sX K — \ (that is, it is not renormalized) and vanishes 
as K — >■ 2~. This vanishing of a is also responsible for 
the SF-MI transition observed for weak lattices in Ref |6j 
at commensurability; K ^ 1 corresponds to a critical 
7c ~ 3.5. We obtain a for intermediate values \ < K < 
2 by direct numerical integration of the flow equations 
given in Ref. [18]. 

To study the system of interacting kinks, we employ 
the statistical density matrix in imaginary time and po- 
sition representation [121 HO] • This is given by 



p{R,R-p)^ n{Bx,eM-S[{xn}]/h) 



(5) 



where {a;„(T)} denotes the set of positions of the parti- 
cles at time t, R ~ {ri, r2, . . . , tat} denotes the set of 
positions of the particles at r = and t = P (see below) 
while Dxj denotes functional integration over Xj] finally. 



S[{Xn}] = I 



dT 



n.m 



V{Xr, 



(6) 

In Eq. the integral runs from zero to /3 and there is an 
ultraviolet cutoff A,- = 27r/AT with At a discretisation 
step size [20]. In this picture, the worldlines of the par- 
ticles Xi{T) correspond to classical strings without over- 
hangs, the ends of which are fixed to 3:^(0) = Xi{T) = ri. 
Note that Eq. ^ and Eq. ^ describe the p appro- 
priate for distinguishable particles; for bosonic parti- 
cles, one symmetrises in the end, so that pb{R, R' ; (3) = 
p{R, PR'; /3), with P labelling the permutation. 
We begin by estimating the temperature dependence 
of the critical incommensurability Qd above which ex- 
change effects become important. The worldlines of the 
particles are of length p in the time-like direction, and 
the "width" of the path in the space-like direction will 
be u; oc ^/Wpja. If the average inter-kink distance, pro- 
portional to Q~^, is larger than this, quantum effects are 
not important; the condition for the statistics to be im- 
portant is therefore QyJWpJa > 1, up to a numerical 



factor. This defines a critical Qd oc yJkBTa/h?. Be- 
low this Qd, the kinks behave like noninteracting bosons; 
above it, they begin to interact, and we expect the ef- 
fects described below to be evident. Furthermore, since 
a vanishes on the lines V = {) and K = 2, Qd also van- 
ishes there. 

Next, we concentrate on the T = or /3 — > oo limit, 
corresponding to infinitely long strings; in this limit, the 
degeneracy condition is always satisfied. We shall employ 
a renormalization group (RG) technique applied directly 
to the density matrix of Eq. Details of this will be 
presented elsewhere [H]; here, we shall only outline our 
conclusions. 

Splitting the fields Xi into slow and fast parts as 
usual [22], it is possible to extend the Wegner-Houghton 
approach [231 121] to the many-body case, obtaining the 
fiow equation for the potential 



d,V = V 



lxd:,V + hAr \0g(l+^dlV 

I \ a At- 



(7) 



Notice that the coarse-graining is done in the r direction; 
thus, information on lengthscales comparable to the kink 
density is still present in the fixed point-potentials. 

The physics of the system is determined by the fixed 
point potentials of Eq. ([t]). For bare (initial) potentials 
that diverge at the origin [22], these may be determined 
numerically; for x 3> h/ K^-a, their behaviour is 



fp 



fi2 A(A - 1) 
2a ' 



(8) 



where we have written the coefficient of x~'^ as h^X{\ — 
l)/2tT in order to make contact with the conventional 
notation in the literature (see below). We are thus 
dealing with a system of bosons interacting via an in- 
verse square power law; this is the celebrated Calogero- 
Sutherland model [10], the ground-state wavefunction 
and low-energy spectrum of which are known. We con- 
centrate here on its ground-state properties, which have 
been studied using numerical techniques [2S]. The au- 
thors of Ref. [21] find (quasi-)long-range off-diagonal or- 
der for < A < 2, while they find (quasi-)long-range 
diagonal order for A > 1. The system is thus in a con- 
densed, SF state for < A < 1, in a SS state, char- 
acterised by the simultaneous presence of diagonal and 
off-diagonal long-range order for 1 < A < 2, and in a 
crystalline, S state characterised by strong diagonal cor- 
relations for 2 < A. Therefore, the phase in which the 
system is for incommensurate densities [Q ^ 0) depends 
on the range in which the A corresponding to the poten- 
tial in Eq. Q hes. 

To map out the phase diagram, we note that local 
analysis of the fixed point ordinary differential equation, 
Eq. ([t]) with the left hand side set to 0, indicates that 
the fixed point potentials, P/p, have the property that 
dVfp/da > (for all x). An increase in a therefore re- 
sults in an increase in A of Eq. ((sj) . In addition, at if = 1 
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and V — (hard-core free bosons), A = 1 26J. Based on 
these two pieces of information, and the behaviour of the 
mass described earher, we propose the phase diagram in 
Fig. [l] for T = 0. Starting from the point K ^ 1, V ^ 0, 
an increase of V causes a rapid increase of a, which cor- 
responds to an increase in A so that A > 1 which corre- 
sponds to SS. As V is further increased, A reaches the 
value A = 2 at y = Vc^m, at which point phase coherence 
is lost, the structure factor displays a sharp peak 
and the system is in the S state. On the other hand, 
starting from any point on the K = 1 line and increasing 
K corresponds to decreasing cr, thus decreasing A from 
its value at K ~ 1. As a result, the line Vc^m curves 
upwards as K increases. Starting from Vcfl{K = 1) = 
and increasing K, A must decrease below 1 so that the 
diagonal order is suppressed; thus, the line 14,0 separat- 
ing the SF from the SS region also curves upwards. As 
-fC 2, or 7 — 3.5, the effective mass vanishes for any 
V; this results in a rapid decrease of A, so that both lines 
curve upwards sharply. 

It is important to note that the presence of the SS 
phase represents an order out of disorder effect: quan- 
tum fluctuations, which at first sight one would expect 
disorder the system, result in a strengthening of the re- 
pulsion which in turn causes the system to order. 

Let us briefly discuss the differences between the 
phases just described in terms of experimentally acces- 
sible quantities. The main distinguishing features of 
these phases are the diagonal and off-diagonal correla- 
tions. Off-diagonal long-range order may be observed 
using time-of- flight measurements, which therefore allow 
us to distinguish the SF and SS phase-coherent phases 
from the S phase; in the latter, phase correlations drop 
quickly and the time-of-flight image is smeared. On the 
other hand, techniques for measuring density variations 
would distinguish between the SS and S phases on one 
hand and the SF phase on the other; single-site address- 
ability is possible [57] , which may alUow to detect density 
oscillations. 

In summary, we have shown that the 
incommensurability-induced MI-SF transition oc- 
curs for arbitrarily small incommensurability. We have 
then studied the system of bosonic quasiparticles which 
appears as soon as commensurability is lost; calculating 
the effective interactions between them, as well as their 
effective mass, and using an RG transformation, we have 
argued that quantum fluctuations drive the interactions 
to acquire an inverse square form (a Calogero-Sutherland 
model). A phase diagram for the current system of 
strongly interacting bosons in a weak ID optical lattice 
is then proposed, which features SF, SS, and S phases. 
The most striking feature of our calculations is the 
appearance of a SS phase, in spite of the local character 
of the interactions in the original model. The periodicity 
of both the SS and the S phases is unrelated to that 



of the underlying lattice, thus providing us with more 
exotic states of matter. 

Although the experimental setup of Ref. [6^ allows, 
in principle, to tune the density and to investigate also 
the commensurate- incommensurate quantum phase tran- 
sition, up to now only the commensurate regime has been 
studied. Our studies indicate that the elusive SS phase 
is within reach by modifying a single parameter in the 
experimental setup of Ref. [6 . We hope that our work 
will trigger further experiments into this fascinating and 
largely unexplored regime. 
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